Spectral Properties of a Quantum Impurity in d-Wave Superconductors 
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The variational approach of Gunnarsson and Schonhammer to the Anderson impurity model is 
generalized to study d-wave superconductors in the presence of dilute spin-1/2 impurities. We show 
that the local moment is screened when the hybridization exceeds a nonzero critical value at which 
the ground state changes from a spin doublet to a spin singlet. The electron spectral functions 
are calculated in both phases. We find that while a Kondo resonance develops above the Fermi 
level in the singlet phase, the spectral function exhibits a low-energy spectral peak below the Fermi 
level in the spin doublet phase. The origin of such a "virtual Kondo resonance" is the existence of 
low-lying collective excitations in the spin-singlet sector. We discuss our results in connection to 
recent spectroscopic experiments in Zn doped high-T c superconductors. 

PACS numbers: 74.25.Jb, 71.27.+a, 72.15Hr, 74.25Ha 
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The physics of quantum impurities in high-T c super- 
conductors has generated an increasing amount of exper- 
imental and theoretical interests. NMR measurements in 
YBCO showed that replacing the Cu-atom in the CuC>2 
plane by a nonmagnetic Zn or Li atom induces a spin-1/2 
magnetic moment on the Cu ions in the vicinity of the 
impurity 1,2 . Atomic resolution scanning tunneling micro- 
scope (STM) experiments on BSCCO found that the lo- 
cal density of states (DOS) near the Zn impurity exhibits 
a sharp resonance peak just below the Fermi energy 3 . 
The origin of the resonance peak has been attributed 
to potential scattering 4,5 , local moment formation 6 , and 
magnetic Kondo scattering 7-10 where it is identified with 
the Kondo resonance arising from the screening of the 
local moment by the electrons in d-wave superconduc- 
tors (dSC). While it remains controversial as to which 
of the above is the dominant mechanism responsible for 
the low-energy conductance peak, and how the findings 
of NMR are reconciled with those of STM, it has become 
important to acquire basic understandings of the spectral 
properties of quantum impurities in dSC. 

The ground state and the impurity spectral function 
in the non- superconducting phase with a vanishing den- 
sity of states (DOS), i.e. the pseudogap phase, have 
been studied in detail following the work of Withoff and 
Fradkin 11,8 . In this paper, we study these properties 
associated with the Kondo physics in the d-wave su- 
perconducting phase. In addition to the ground state 
properties, we determine the momentum dependent, con- 
duction electron spectral function A{k,uS) in dSC. This 
is an important quantity since it is directly measured 
by angle-resolved photoemission spectroscopy (ARPES). 
To this end, we describe the localized impurity by the 
infinite-U Anderson model and generalize the variational 
large-N approach of Gunnarsson and Schonhammer 12,13 
for rare earth and actinide heavy fermion materials to 
the case of dSC. This is a reliable, physically transpar- 
ent many body approach that allows an essentially ana- 
lytical calculation of A(k,Lu). Specifically, starting from 
the singlet BCS state with the impurity site empty (sin- 



glet) or singly occupied (doublet), we construct the vari- 
ational wave functions of states in both the spin singlet 
and doublet sectors from a basis of states generated by 
the hybridization between the impurity and the conduc- 
tion electrons to next to leading order in 1/N. We find 
that there exists a critical value for the hybridization 
that separates a spin doublet ground state from a singlet 
ground state where the local moment is fully screened. 
Using the variational wave functions, we calculate the T- 
matrix, the self-energy, and the spectral function of the 
conduction electrons in both the singlet phase and the 
local moment phase. We find that while a Kondo res- 
onance naturally develops above the Fermi level in the 
singlet phase, A{k,ui) exhibits a spectral peak of a simi- 
lar strength below the Fermi energy in the local moment 
phase. The origin of this virtual Kondo resonance can be 
traced to the existence of low-lying collective spin singlet 
excitations above the doublet ground state in the local 
moment phase. The qualitatively different behaviors of 
the spectral function make it possible to distinguish be- 
tween the two ground states by STM and ARPES. 

The single impurity Anderson model in a dSC can be 
written as H = H c + H d + H v , 



ka k 

H d = e d J2 dUe + Ud\d ] d\d i 

a 

H v =Vj2(cLd a + h.c). 



(1) 
(2) 
(3) 



Here c ka creates an electron with a dispersion et — 
cos k y ) — /i in a dSC described by the gap 
= iA (cosfc a; — cos k y ); S a creates a local- 



— 2t(cosk x -\ 
function 

ized electron with energy e d at the impurity site, and U 
is the on-site Coulomb repulsion. We take the infinite-U 
limit such that the impurity site is either empty or singly 
occupied. The hybridization between the impurity and 
the conduction electrons is described by the hybridization 
constant V in H v . For simplicity, only the case of N = 2 



1 



is written explicitly in Eq. (1-3). The generalization to 
the case of N orbital degeneracy with NV 2 — 0(1) is 
straightforward 12-14 . 

As in the BCS theory, H c can be diagonalized by 



a Bogoliubov transformation: c 



and cl kl = 



u k a 



t 



u k a 



+ v k a kh where u\ = 1 - v\ = 

(1 + e k /E k )/2 with E k = ^el + Al. H c and H v can 
be expressed in terms of the quasiparticles, H c = 

T,ka E ka{ a a k a, and 
k 

+ ( u k^_ kl +v k a k t)d l +h.c]. (4) 
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FIG. 1. Schematics of the variational states. Singlet 
(first row) and doublet (second row) states to leading (a,b, 
and e) and next to leading order (c,d,f, and g) in 1/N. 

We now generalize the variational method of Gunnars- 
son and Schonhammer to the case of superconductors 15 . 
Since the BCS ground state \BCS) with the impurity 
level empty is a singlet (Fig. la), we can construct other 
states in the singlet sector by acting on \BCS) with 
the Hamiltonian H . Operating once with the hybridiza- 
tion term generates the type of states \edoE k <j) shown in 
Fig. lb in which a quasiparticle excitation with spin a is 
created from the dSC condensate and the impurity level 
is singly occupied by an electron with the opposite spin 
a. Taking the linear combination of all such states gives 
the trial wavefunction for the singlet state: 

\S) = A s [\BCS) + -L £ Bj^ \wE kl <r) 

fcl<7 



+ 



-y^ C k lk2 \E kl (rE k2 a) 



E 



y - 

N 



D s klk2k3 \e d a'E kl aE k2 aE k3 a')]. (5) 



The last two terms in Eq. (5) contain singlet states 
generated to the next to leading order in the 1/N ex- 
pansion shown in Figs, lc and Id. They correspond 
to states with two quasiparticle excitations and an un- 
occupied impurity level and three quasiparticles plus a 
singly occupied impurity level respectively. The param- 
eters A s to D s are determined by minimizing the energy 
E s = (S\H\S)/(S\S). We find 

B{ = VNVv k [E s -e d -E k - ^(E s - E k )}-\ 



C s klk2 = Vu k2 B s ki [E s - E 12 - £ (£ s - e d - E 12 )}~\ 
D kk 2 k 3 = ^NVv k3 C s klk2 [Es -e d - E 123 }-\ (6) 
where E 12 = E kl + E k2 , E 123 = E kl + E k2 + E k3 , and 



NV 2 v\ 



Es-ed-Ek-S^Es-EkY 



(7) 
(8) 



Z (u) = NV 2 J2vl(u-E k )-\ 

k 

Si(w) = V 2 J2 U > - Ek - So(w -ea- EkT 1 . (9) 
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FIG. 2. Phase diagram of the infinite-U Anderson 
model in a dSC to leading (dashed line) and next to lead- 
ing order (solid line) in 1/N. 

The variational state in the spin doublet sector can be 
constructed starting from the BCS state plus a singly oc- 
cupied impurity level (Fig. le) in the large-N limit. To 
next to leading order in 1/N, two types of states are gen- 
erated: one corresponds to a quasiparticle excitation and 
an empty impurity level (Fig. If) and the other to two 
quasiparticles out of the condensate and a singly occu- 
pied impurity level shown in Fig. lg. Denoting the net 
moment in the doublet state by r, we have 

\Dt) = A d [\BCSe d T)+J2B d k \E k T) 

k 

+ 7^ E C d klk2 \e d aE kl rE k2 a)]. (10) 

As in the singlet case, the variational parameters A d , B d , 
C d , and the total energy are determined by minimizing 
E d = (Dt\H\Dt)/(Dt\Dt) in the doublet state: 

Ek = Vu k [E d — E k — ^o(E d — E k — e d )] x , 
Ct k2 = VNVv k2 B kl (E d - E kl - E k2 - e d )-\ 



V 2 u\ 



Ed - E k - S (E d — Ek — ed) 



(11) 
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The ground state and the phase diagram are obtained 
by comparing E s and E d derived in Eqs. (7) and (11) as 
a function of the impurity level e d and the hybridization 
V. We use the d-wave gap A as the energy unit, set t 
to A , and integrate over k in the continuum limit. We 
find that, due to the linearly vanishing DOS in the dSC, 
the impurity spin can be fully screened by the Bogoli- 
ubov quasiparticles only if the hybridization V is greater 
than a nonzero value. Similar to the pscudogap Kondo 
impurity model 8 ' 11 , both the spin singlet and the doublet 
local moment phases exist for an Anderson impurity in a 
dSC and are separated by a line of phase transitions in 
the phase diagram shown in Fig. 2. 

One of the advantages of the variational wave func- 



tion approach is that it allows a reliable and transpar- 
ent many-body calculation of the spectral function 12 , 
A(k,uj) = (— l/27r)ImG c (fc, uj), where G c (fc,w) is the 
impurity averaged, retarded normal Green's function 
for the conduction electrons. The local Green's func- 
tion G c (i,j,T,t) = -i9(t)(^ \{c ia (t),c f j!T (0)}\(l)o) where 
\4>o) — \S), \Dt) represent the fully interacting singlet 
and doublet (with polarization r) ground states. Both 
the local electron Green's function at the impurity site 
and the impurity averaged conduction electron Green's 
function in momentum space can be written in terms of 
the T-matrix, which can be obtained from the variational 
wave functions. Using the Nambu representation, the lo- 
cal Green's function can be written as: 



G c (r, r', uj) = G°(r - r', uj) + G°(r, w)f (u,)G°(-r', w) f(w) = \V\ 2 G d (co), 



(12) 



where Gd(u>) is the local electron Green's function at the impurity site. The impurity averaged Green's function is 
given by 



G c (k,u>) 



(13) 



where p lmp is the density of impurities. Using the variational wave functions, G d (uj) can be derived in closed form to 
first order in 1/N. For the singlet ground state we obtain, 



G s d ^w + iS) = (S\di 



uj + id - H + E s 



d\\S) + (S\d\ 



1 



w-e d - T, d {tJ + iS + E s ) + E s + iS 
1 

u> + e d + H d (-w -i6 + E s ) - E s + iS 




d T |S) 



u k V 



N K uj + E k - E s + iS 



lu + E k - E s + iS 



(14) 



where Sd(w + iS) — J2k U V^ /( w ~ + iS). In the local moment phase, the spin SU (2) symmetry is broken and the 
STM measures the spin averaged tunneling density of states. Since the total Hamiltonian is invariant when the spin 
up and down electrons are interchanged, averaging over the conduction electron spin is equivalent to averaging over 
the polarization of the local moment ground state. We thus obtain 



G£ (w + iS) = \g% (w + iS) + \G% (w + iS) = 



r=1,l 

A 2 r 



1 



uj + iS - H + E D 
1 



d\\Dr) + (Dr\d\ 



1 



uj + iS + H - E D 

v k V 



d T |Dr) 



uj + Eo(— w — i5 - ed + Ed) - Ed +i5 
1 



\ k 



uj - E (w + iS - e d + E D ) + E D +iS 



uj - ed - Ek + Ed + id 

i V — B d2 - 

^ /v k ,,, - ej - F„_ 



u-e d - E k + E D +iS 



(15) 



We next present results for the spectral function of 
the conduction electrons. The impurity level is fixed at 
e<j = — 1.4A . We choose V = 1.6A for the doublet local 
moment phase at a corresponding valence (rid) = 0.93, 
and V — 1.8A for the singlet state with (n d ) = 0.58. 



We will show that along nodal directions of the d-wave 
gap, where SC coherence is absent (it 2 . = 1, v\ = for 
\k\ < kf), the spectral function exhibits qualitatively dif- 
ferent behaviors in the two different phases. In Figs. 3a 
and 3b, the solid lines depict the low-energy spectra along 
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the nodes at k x = k y = 0.63/c/. In the singlet phase 
(Fig. 3b) , a resonance peak appears above the Fermi level, 
which is the expected Kondo resonance. Physically, this 
can be understood by considering an intermediate state 
in the infinite-U limit with one less electron than the sin- 
glet ground state. An added electron at the resonance 
energy creates a quasiparticle excitation and forms a sin- 
glet state that has a large overlap with the Kondo sin- 
glet state. The intrinsic width of the Kondo resonance 
is captured by the present variational wave function ap- 
proach. In contrast, the spectrum in the local moment 
phase (Fig. 3a) shows, remarkably, a resonance peak be- 
low the Fermi level. This peak does not have the usual 
meaning of a Kondo resonance. Instead, it can be under- 
stood as a virtual Kondo resonance. Although the ground 
state is a doublet and the local moment is unscreened, 
the collective singlet states of the type given by Eq. (5) 
exist as low energy excitations above the doublet ground 
state. When an electron is removed from the SC con- 
densate, quasiparticle excitations with the opposite spin 
can be created to screen the local moment and occupy 
the virtual spin singlet state. Note that, in general, de- 
pending on the sign of the particle-hole asymmetry, the 
Kondo resonance in the singlet phase can be either above 
(the present case for the infinite-U Anderson model) or 
below the Fermi level. The virtual Kondo resonace in the 
local moment phase is located on the opposite side of the 
Fermi level compared to the Kondo resonance. 
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FIG. 3. Momentum resolved spectral functions in the 
local moment (a) and the singlet (b) phases along the 
nodal (solid lines) and anti-nodal (dashed lines) direc- 
tions. 



The behavior of the spectra along anti-nodal directions 
of the d-wave gap is shown in dashed lines in Figs. 3a and 
3b at k x = kf, k y — 0. Due to the strong SC coherence 
{u\,vj, 7^ 0), the Kondo and the virtual Kondo reso- 
nances are reflected on the hole and the particle sides 
such that the resulting spectra show nearly symmetric 
resonance peaks. These results suggest that both ARPES 
and STM experiments STM can be used to distinguish 
between these two ground states. 

Finally, we turn to the local density of states (LDOS) 



as measured by STM experiments, 

p(r,T,oj) = ImGJk'a, ka, r, oj)e 

k,k' ,a 



-i(k — k')r 



In Fig. 4, the LDOS at the impurity site is shown for 
both the singlet and the local moment phases. Due to 
the high symmetry of this tunneling point, the LDOS 
spectrum is dominated by the contribution from the low 
energy quasiparticle excitations along the nodal direc- 
tions. As a result, the tunneling spectrum closely resem- 
bles that of A(k, uj) along nodal directions: a single reso- 
nance peak above the Fermi level in the singlet phase and 
below the Fermi level in the local moment phase. The 
STM data 3 interpreted this way would suggest that the 
Zn-doped BSCCO surface is in the local moment phase. 
The LDOS away from the impurity is shown in the in- 
sets of Fig. 4. The conductance spectra in both phases 
still exhibit a single peak along nodal directions, whereas 
along antinodal directions two peaks appear due to the 
SC coherence. The amplitude of the peaks is significantly 
reduced than at the impurity site. 




FIG. 4. The LDOS spectrum at the impurity site in 
the singlet (dashed line) and the local moment (solid line) 
phases. The insets show the LDOS away from the impu- 
rity along nodal (dashed lines) and anti-nodal (solid lines) 
directions in local moment (b) and singlet (a) phases. 

To conclude, we have studied the spectral properties 
of electrons in dSC coupled to localized Anderson im- 
purities in both the singlet and the local moment phases 
using the variational wave function approach. Our results 
show that Kondo screening is not a prerequisite for the 
emergence of low energy resonance peaks in the spectral 
function. Even if the impurity spin is not screened, the 
virtual Kondo resonance due to collective spin singlet ex- 
citations would still lead to a sharp resonance peak, but 
located on the opposite side of the Fermi level compared 
to the Kondo resonance in the singlet phase. It is the 
evolution of the resonance peak positions that reveals the 
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screening properties of the local moment in the ground 
states, which can in principle be measured by ARPES 
and STM experiments. 
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